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ABSTRACT 



Several papers have recently highlighted the possibility of measuring redshift space 
distortions from angular auto-correlations of galaxies in photometric redshift bins. In 
this work we extend this idea to include as observables the cross-correlations between 
redshift bins, as an additional way of measuring radial information. We show that 
this extra information allows to reduce the recovered error in the growth rate index 
7 by a factor of ~ 2. Although the final error in 7 depends on the bias and the 
mean photometric accuracy of the galaxy sample, the improvement from adding cross- 
correlations is robust in different settings. Another factor of 2 — 3 improvement in the 
determination of 7 can be achieved by considering two galaxy populations over the 
same photometric sky area but with different biases. This additional gain is shown to 
be much larger than the one from the same populations when observed over different 
areas of the sky (with twice the combined area) . The total improvement of ~ 5 implies 
that a photometric survey such as DES should be able to recover 7 at the 5—10% from 
the angular clustering in linear scales of two different tracers. It can also constrain the 
evolution of f(z) x a%(z) in few bins beyond z ~ 0.8 — 0.9 at the 10—15% level per-bin, 
compatible with recent constrains from lower-z spectroscopic surveys. We also show 
how further improvement can be achieved by reducing the photometric redshift error. 

Key words: galaxy clustering; redshift-space distortions; angular correlations; pho- 
tometric redshift surveys 



Our understanding of the local Universe and the way it 
evolved from small perturbations has been reshaped over the 
past decades with the successful completion of vast observa- 
tional campaigns for CMB fluctuations, large scale structure 
and SNIa distances. Yet several still open issues arose from 
these studies, the most important of which is probably the 
late-time accelerated expansion of the Universe. 



1 INTRODUCTION 





Redshift space distortions (RSD) (Kaiser 1987; Hamil- 
ton 1998) can be used to understand the (linear) growth of 
structures, which provides a direct path to study the origin 
of cosmic acceleration. On large scales, RSD arises from the 
coherent velocities of galaxies and reveals how perturbations 
grow in time. Typically this method requires measuring of 
galaxy clustering in 3 dimensions (3D) in order to sample di- 
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rections parallel and transverse to the line-of-sight where the 
effect is maximized or cancels out completely (see e.g. Oku- 
mura 2008; Guzzo et al. 2009; Cabre & Gaztanaga 2009; 
Blake et al. 2011; Reid et al. 2012; Kazin et al. 2013 and 
references therein). 

Over the past few years it has been however shown that 
the effect of RSD is also present, albeit with a smaller contri- 
bution, in the angular (2D) clustering of photometric galaxy 
samples if they are selected in photometric redshift bins, see 
for instance Nock, Percival,& Ross 2010; Crocce, Cabre, & 
Gaztanaga 2011; Ross et al. 2011. This concrete idea has 
been already applied to data using a sample of photometric 
Luminous Red Galaxy (LRG, see Padmanabhan et al. 2007; 
Thomas, Abdalla & Lahav 2010; Crocce et al. 2011). 

Yet all the previous studies focused on the angular clus- 
tering from a set of measurements of auto-correlation in one 
or more redshift bins. The goal of this paper is, on the one 
hand, to extend these analysis to include also the cross- 
correlations between redshift bins in order to account for 
some radial information. This is motivated by the recent 
findings of Asorey et al. 2012 who show how a tomographic 
(2D) study involving auto and cross correlations can yield 
the same constrains on cosmological parameters as a full 
spatial (3D) study. 

On the other hand we will also investigate the improve- 
ments brought by considering two different populations (and 
their cross correlations) in the likelihood analysis for the 
growth rate. This is motivated by the fact that for the spec- 
troscopic analysis, the combination of different samples trac- 
ing the same underlying matter fluctuations can be used 
to decrease sampling variance and improve considerably the 
constrains in growth of structure (McDonald & Seljak 2009; 
White, Song & Percival 2009; Gil-Marin et al. 2010). 

This paper is organized as follows. In Sec. [2] we lay out 
the methodology, including the analytical tools, the defini- 
tion of the different samples and surveys and the likelihood 
analysis. In Sec. [3]we present our result, and in Sec.[4]our 
conclusions. 



2 METHODOLOGY 

Our goal is to study the effect of RSD in angular clustering, 
especially it usefulness to derive constrains on the growth 
of structure at large scales. We study angular clustering us- 
ing auto- and cross- correlations between redshift bins. The 
inclusion of cross correlations between different radial shells 
allow us to include the radial modes that account for scales 
comparable to the bin separation. On the other hand, the 
angular spectra of each redshift shell includes information 
mainly from transverse modes. 

With the idea of a potential sample variance mitigation 
in the analysis, we also consider the correlation between the 
angular clustering of different tracers of matter, considering 
them either independent (i.e. each tracer in a different patch 
of the sky) or correlated (same sky). 

Throughout this paper we use CAMB_source (Lewis 

B camb. info/sources 



Number of bins N z 


Az/(l + z) 


4 


0.15 


6 


0.1 


8 


0.08 


12 


0.05 


19 


0.03 



Table 1. The different redshift bin configurations considered in 
our paper, within a photometric redshift range of 0.4 < z < 1.4. 
We show the total number of bins and their redshift width Az 
(which evolves with redshift in the same manner as the photo-z) . 

et al. 2000; Lewis et al. 2007; Challinor & Lewis 2011), 
including cross correlations between radial bins and the cor- 
relations between different populations. Let us note that we 
use the exact Ct computation in CAMELsources, because in 
angular clustering the imprint of redshift distortions affect 
mainly the largest scales, which are not included when us- 
ing the Limber approximation (Limber 1954; Lo Verde & Af- 
shordi 2008; Crocce, Cabre, & Gaztanaga 2011). Moreover 
the Limber approximations does not account for clustering 
in adjacent redshift bins. 

2.1 Fiducial survey and galaxy samples 

We start by describing the fiducial photometric survey that 
we assume in our analysis (characterized by a redshift range 
and a survey area) and the different galaxy samples consid- 
ered within that volume (characterized by the bias b, the 
accuracy of photometric redshift estimates a z and their red- 
shift distribution). 

Our fiducial survey is similar to the full DES, with an 
area coverage of one octant of the sky (i.e., f s k y = 1/8) and 
a redshift range 0.4 < z < 1.4. We characterize the redshift 
distribution of galaxies within this survey by 

dzdtl 9al \0.5J (> 

where Ng al is a normalization related to the total number 
of galaxies of each population sample, denoted by a. We 
typically consider two types of sample populations, one with 
bias b = 1 and a z = 0.05(1 + z) (Popl) and another with 

6 = 2 and a z — 0.03(1+2) (Pop2). For simplicity we consider 
the same redshift distribution for all samples with a fiducial 
comoving number density, unless otherwise stated, of n(z — 
0.9) = 0.023/i 3 Mpc~ 3 . This value corresponds to a total of 
~ 3 x 10 s galaxies within the survey redshift rang^] 

In Table [l] we show the different redshift binning 
schemes in which we divide our survey prior to study the 
clustering either with the auto-correlations or with the 2D 
tomography that also includes the cross-correlations be- 
tween bins. Note that we consider consecutive bins with 
an evolving bin width with redshift, i.e. Az oc (1 + z), to 
match the photometric uncertainty which also assumes a 
linear evolution with redshift. 

7 Note that this matches the nominal 3 X 10 8 galaxies expected 
to be targeted above the magnitude limit of DES (i < 24) 
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2.2 Angular power spectrum 

In our analysis we study angular clustering using the angular 
power spectrum of the projected overdensities in the space of 
spherical harmonics. The auto-correlation power spectrum 
at redshift bin i, for a single population, is given by: 



( -," = 2 / dk k 2 p (k) (*i(k) + 



*j(fc)= / dz <t>i(z)b(z)D(z)j e (kr(z)) 



where 



is the kernel function in real space and 

2l 2 + 21 - 1 



(2) 



(3) 



dz 4> l (z)f(z)D{z) 
l{£-l) 



{21- 1)(2£+1) 

(e + i)(e + 2) 



(2£ + 3)(2£- 1) 



ji(kr) 



2£ + l)(2£ + 3) 



je+2{kr) 



(4) 



should be added to if we also include the linear Kaiser 
effect (Fisher, Scharf & Lahav 1994; Padmanabhan et al. 
2007). In Eqs. | 3|4 | b(z) is the bias (assumed linear and 
deterministic), D(z) is the linear growth factor and f(z) = 
din D/d\na is the growth rate. Photo-z effects are included 
through the radial selection function 4>{z), see below. 

For the case of 1 population, there are N z auto- 
correlation spectra, one per radial bin. Then, we add to our 
observables the N Z (N Z — 1)/2 cross-correlations between dif- 
ferent redshift bins. These are given by 

C i = lf dk fc2jP W + *< *"(*)) (*<(*) + r (*) 

(5) 

Therefore, we are considering N Z (N Z + l)/2 observable an- 
gular power spectra when reconstructing clustering informa- 
tion from tomography using N z bins, for a single tracer. 

If we combine the analysis of two tracers, a and /?, the 
angular power spectrum is given by 



= - J dk k 2 P(k) (^l a (k) 

+ (*^(fe)+*f r (fe)), 



(6) 



where ^\ and characterize each galaxy sample through 
the radial selection function <pi(z) and the bias b(z) in ex- 
pressions and @ . We use the general notation where 
C^""'' 3 is the correlation between redshift bin i of population 
a with redshift bin j of population (3. By definition, 



C 



*3fS 



- U £ 

± for a^fci^j 



(7) 
(8) 



Then the total number of observables is 2N Z (2N Z + l)/2 if 
we consider the same redshift bins configuration for both 
populations, in the case in which both are correlated. 



2.2.1 Radial selection functions 



The radial selection functions (f>i in Eqs. ( 2|5|6 1 encode the 
probability to include a galaxy in the given redshift bin. 
Therefore, they are the product of the corresponding galaxy 
redshift distribution and a window function that depends on 
selection characteristics (e.g binning strategy), 



dz 



(9) 



where dN a /dz is given by Eq. Q. We include the fact that 
we are working with photo-z by using the following window 
function: 



Wi(z) 



dz p P(z\z p )W[ h (z p ) 



(10) 



where z p is the photometric redshift and P(z\z p ) is the prob- 
ability of the true redshift to be z if the photometric estimate 
is z p . For our work we assume a top-hat selection W[ h (z p ) 
in photometric redshift and that P(z\z p ) is Gaussian with 
standard deviation a z . This leads to, 



dN a 
dz 



erf 



erf 



V2af 



(11) 



where z p , m i„ and z Pl max are the (photometric) limits of each 
redshift bin considered and a z is the photometric redshift 
error of the given population a at the corresponding redshift. 



2.2.2 Covariance matrix of angular power spectra 

We assume that the overdensity field is given by a Gaussian 
distribution and therefore, the covariance between correla- 

is given by 



tion C^"' 7 ' 3 and correlation G\ 



Cov^ ( 



c: 



^0 



N(l) 



(12) 

where N(l) = (21 + l)A£f a k y is the number of transverse 
modes at a given I provided with a bin width A£. We set 
A£ = 2/fsky, the typically chosen value to make Cov block- 
diagonal (Cabre et al. 2007; Crocce, Cabre, & Gaztanaga 
2011). In this case, bins in £ are not correlated between them. 

Therefore, for each I bin, we define a matrix with 
2N Z (2N Z + l)/2 rows, where N z is the number of red- 
shift bins, taking into account the covariances and cross- 
covariances of auto and cross-correlations between each pop- 
ulation and among them. In order to include observational 
noise we add to the auto-correlations of each population in 



Eq. ( 12 1 a shot noise term 



C 



= c 



An 



(13) 



that depends on the number of galaxies per unit solid an- 
gle included in each radial bin. We define the \ assuming 
the observed power spectrum Cf" 1 corresponds to our fidu- 
cial cosmological model discussed in Sec. < |2.3[ ), while we call 
C{ nod the one corresponding to the cosmology being tested, 



X 



- E [pt 



Cmod 
I 



Cov," 1 C, 



C? od ). (14) 
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Figure 1. The gain from adding redshift-bins cross-correlations. Dashed lines show the expected 1-a constrains in 7 from the combined 
analysis of angular auto-correlation in photo-z bins spanning 0.4 < z < 1.4, as a function of the bin width Az/(1 + z) (see Table[l]for 
the corresponding total number of bins). Different colors correspond to different populations with bias and u z as labeled. Solid lines 
show, for each population, the same study but also including all the cross-correlations between bins (and their complete covariance). For 
optimal bin widths Az < cr z the gain from including cross-correlations is ~ 2 or better. 



2.3 Cosmological model and growth history 

Throughout the analyses, we assume the underlying cosmo- 
logical model to be a flat ACDM with cosmological param- 
eters w = -1, h = 0.7, n s = 0.95, Q, m = 0.25, Q b = 0.045 
and erg = 0.8. These parameters specify the cosmic history 
as well as the linear spectrum of fluctuations Pq. In turn, 
the growth rate can be well approximated by, 



f(z) = n m ( z y 



(15) 



and 7 = 0.545 for ACDM. Consistently with this we obtain 
the growth history as 



D[z) = exp 



/(*) 
1 + z' 



dz 



(16) 



(where D is normalized to unity today) . The parameter 7 is 
usually employed as an effective way of characterizing mod- 
ified gravity models that share the same cosmic history as 
GR but different growth history (Linder 2005). Our fiducial 
model assumes the GR value 7 = 0.545. In order to fore- 
cast the constrains on 7 we consider it as a free parameter 
independent of redshift. 

With these ingredients, we do a mock likelihood sam- 
pling in which we assume that the theoretical values for the 
correlations at the fiducial value of the parameters corre- 
sponds to the best fit position. The likelihood is based on the 



X given in ( 14 1. In our case, we keep fixed all the parameters 



and only allow 7 to vary, and then we estimate 68% confi- 
dence limits of it. In the case in which we show constrains on 



fas, we vary this quantity (that now depends on redshift, 
thus the number of fitting parameters is a function of the bin 
configuration), fixing the rest of parameters. The maximum 
t considered in the analysis is £ max = r(zsurvey)kmax ~ 220 
for kmax = 0.1 hMpc -1 , while for the largest scales we set 
Imin = 2. We had to adapt CAMB_sources in order to con- 
strain 7 or fag using the same technique described in the 
Appendix A of Asorey et al. 2012. 



3 RESULTS 

In this section we discuss the constrains on the growth index, 



7 defined in Eq. ( 15 I as obtained for the different redshift 



bin configurations of Table [T] First of all, we study how well 
we can determine 7 using different single galaxy popula- 
tions but including as observables also the cross correlation 
between bins (for a given single population). We also study 
how the constrains depend on the bias and in the photomet- 
ric redshift accuracy of the different galaxy samples. Then, 
we study the precision achievable when one combines dif- 
ferent tracers in the analysis and how this depends on bias, 
photo-z and in particular, the shot-noise level of the sample. 

Lastly we discuss the constrains that we obtain when 
looking into the more standard f(z)a&(z) as a function of 
redshift, and consider auto and cross-correlations of one or 
two galaxy samples. 
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Figure 2. Dependence on photo-z (top panel) and bias (bottom 
panel) for a one-population constrains in 7, as a function of bin 
width (same as in Fig. [TJ. The panels show that lower b and/or 
lower a z yields better constrains in 7. This is hence a competing 
interplay because lower b would correspond to a fainter sample 
with typically worse photometric errors. The dependence on bias 
seems however slightly stronger. 



3.1 Redshift-space distortions with a single 
photometric population 

Let us first consider the constrains on the growth index us- 
ing single photometric populations. Figure [l] shows the \-a 
errors expected on 7 from a combined analysis of all the 
consecutive photometric redshift bins in the redshift range 
0.4 < z < 1.4 as a function of the bin width (i.e. each of the 
configurations detailed in Table 1 J] 

In red we show the constrains on 7 corresponding to 
an LRG-type sample, with bias 6 = 2 and a photometric 
redshift c z /(l + z) = 0.03 (Pop2). Blue lines correspond to 
an unbiased population with cr z /(l + z) = 0.05 (Popl). 

Dashed lines correspond to the case in which we only 
use the auto-correlations in each redshift bin while solid 
lines corresponds to the full 2D analysis that includes all 
the cross-correlations in our vector of observables. Recall 
than in the first case the cross-correlations are included in 
the covariance matrix of the auto-correlations (but not as 



8 Note that different redshift bins can be strongly correlated de- 
pending on bin width and photo-z. We do include this covariance. 



observables). We see that constrains from a full 2D analy- 
sis, including auto and cross-correlations are a factor ~ 2 or 
more better than those from using only auto-correlations. 

From Fig 1 it is clear that in all cases the bin configura- 
tion can be optimized, with the best results obtained when 
Az ~ a z . In addition, there is a competing effect between 
<t z and bias. For broad bins (Az ^> a z ) the photo-z of the 
populations is masked in the projection and the bias domi- 
nates the 7 constrains. Smaller bias gives more relevance to 
RSD and better 7 constrains. As one decreases the bin width 
the population with better photo-z (typically the brighter, 
with higher bias), denoted Pop2, allows a more detailed ac- 
count of radial modes improving the derived errors on 7 
more rapidly than Popl until they become slightly better. 
This optimization is possible until one eventually reaches bin 
sizes comparable to the corresponding photo-z (what sets an 
"effective" width) and the constrains flatten out. 

In Fig. [2] we study in more detail the dependence of con- 
strains with respect to galaxy bias b and photo-z accuracy. 
In the top panel of Fig. [2] we show standard deviation of the 
growth index, A7, fixing the sample bias to 6 — 1 and al- 
lowing two values for photo-z accuracy. Red line represents 
a sample in which u z /(l + z) = 0.05 while blue line has an 
error of cr z /(l + z) = 0.03. In both cases the constrain flat- 
tens once Az ~ a z and the optimal error improves roughly 
linearly with a z . The dependence on the linear galaxy bias, 
6, is shown in the bottom panel of Fig. [2] (for fixed a z ). We 
see that the constrains degrade almost linearly with increas- 
ing bias (see also Ross et al. 2011). As discussed before, this 
is because the lower the bias the larger the relative impact 
of RSD, which results in better constraints on 7. 

In summary we have shown that using the whole 2D 
tomography (auto+cross correlations) allows considerable 
more precise measurements of 7, a factor of 2 or better once 
the bin width is optimal for the given sample. Hence in what 
follows we concentrate in full tomographic analysis. 

3.2 Redshift-space distortions with 2 photometric 
populations 

We now turn to an analysis combining two galaxy popu- 
lations as two different tracers of matter. In Figure [3] we 
compare the constrains from single tracers with respect to 
the combination of both. As before the populations used in 
the comparison correspond to b — 1 and cr z /(l + z) = 0.05 
(Pop 1) and a population with 6 = 2 and <J Z /(1 + z) = 0.03 
(Pop 2). Their respective constrains in 7 are the dashed red 
and blue lines (same as solid lines in Fig. [I]). 

If we combine both tracers and their cross-correlation 
in the same analysis we obtain the constrains given by black 
solid line, notably a factor of 2 — 3 better compared to the 
optimal single population configuration. 

In order to understand how much of this gain is due to 
"sample variance cancellation", in analogy to the idea put 
forward in McDonald & Seljak 2009, we also considered 
combining the two samples assuming they are located in 
different parts of the sky (and hence un-correlated) . We call 
this case Popl+Pop2 in Fig. [3] (solid green line). In such 
scenario the total volume sampled is the sum of the volumes 
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Figure 3. The gain from combining galaxy populations: Comparison of the 68% standard deviations in the growth index from single 
population analysis (dashed lines) with respect to the combined analysis of these two populations over the same field (black solid), using 
all the angular auto and cross-correlations. Remarkably the combination yields errors at least 2 times better than any of the single 
population cases. The solid green line corresponds to the combination of the two samples assuming they are independent (i.e. from 
different parts of the sky). As shown, the combination of correlated populations (same sky) yield stronger constrains than any other case. 



sampled by each population (in our case, two times the full 
volume of DES). This explains the gain with respect to the 
single population analysis. Nonetheless, the "same sky" case 
Popl x Pop2 (where cosmic variance is sampled twice) still 
yields better constrains, a factor of ~ 1.5 — 2, even though 
the area has not increased w.r.t. Popl or Pop2 alone. 

In all, the total gain of a full 2D study with two popula- 
tions (including all auto and cross correlations in the range 
0.4 < z < 1.4) w.r.t. the more standard analysis with a 
single population and only the auto-correlations in redshift 
bins (dashed lines of Fig. [l| can reach a factor of ~ 5. 

As a next step we show how the combined analysis of 
two tracers depends on the relative difference on the bias 
and photo-z errors of the populations. In Fig. [4] we keep 
Popl fix (with b = 1 and a z /(l + z) — 0.05) and we vary the 
bias of Pop2 from 6 = 2 (LRG type bias) to 6 = 3 (galaxy 
clustering like). We keep a z / {1 + z) = 0.03 fixed for Pop2. As 
expected, increasing the bias difference between the samples 
improves the constrains on 7 in a roughly linear way. 

If we now have an unbiased tracer and a highly biased 
one with b = 3, while both tracers have the same <r z /(l + 
z) — 0.03 we obtain constrains given by the black line in 
Fig. [5] Those constrains are better than the case in which 
the unbiased galaxies photo-z is worse, u z /(l + z) = 0.05 
(given by the dashed blue line) . Therefore, if we determine 
photometric redshifts of the unbiased galaxies with higher 
accuracy we will be able to measure the growth rate with 
higher precision. 

One caveat so far is that we have always assumed that 
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Figure 4. Dependence on bias. Increasing the bias difference be- 
tween the samples improves the constrains on 7. The solid black 
line corresponds to the combination Popl X Pop2 of a highly bi- 
ased sample such as LRGs (Pop2) with an unbiased one (Popl), 
while the blue dashed to cluster-like bias tracer as Population 2. 



biases are perfectly known (bias fixed). Hence, in the top 
panel of Fig. [6] we show how the constrains on 7 change if 
we instead consider them as free parameters and marginalize 
over. We see that the difference is very small, in particular 
once the bin configuration is optimal. The reason for this is 
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Figure 5. Dependence on photometric redshift error. Similar 
to Fig. [4] but now changing the photo-z of the unbiased sample 
(Popl) for a fixed 2nd population. The error on 7 depend roughly 
linear with <r z /(l + z) for optimal bin widths. 



clear from the bottom panel that shows the relative error 
obtained for the bias of each sample in the bias free case. 
Because the bias is so well determined (sub-percent) the 
marginalization over them does not impact the error on 7. 

One further concern in our results is that we have as- 
sumed a perfect knowledge of the galaxy redshift distribu- 
tions for both samples. In a more realistic scenario the dis- 
tribution of photometric errors will be known with some 
uncertainty. In order to study the impact of this potential 
unknown we repeated the PoplxPop2 case (for a bin config- 
uration Az/(1 + z) — 0.1) this time marginalizing over the 
value of (Ti z and c?2z when determininig the constrain in 7 
(instead of fixing their values to 0.05(1 + 2:) and 0.03(1 + 2)). 
We find that the resulting A7 only increases by ~ 10% or 
less with respect to the error shown by the black solid line 
in Fig. [3] Such a small change is because the values of o~\ z 
and o"2 2 are very well constrained after the marginalization, 
similar to what happens with the bias free case above. 

An even more thoughtfull analysis of this issue would 
allow for independent errors in the redshift distributions at 
each redshift bin (with some priors) rather than a global 
change of the mean photometric error. And possibly also a 
marginalization over cosmological parameter space. However 
our results above indicate that this should not have a major 
impact in our conclusions. Hence we leave this study for 
future work. 



3.2.1 The impact of shot-noise 

One strong limitation when it comes to implementing the 
"multiple tracers" technique in real spectroscopic data is 
the need to have all the galaxy samples well above the shot- 
noise limit (at the same time as having the largest possible 
bias difference), see for instance Gil-Mann et al. 2010. This 
is cumbersome because spectroscopic data is typically sam- 
pled at a rate only slightly above the shot-noise (to maximize 
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Figure 6. Bias free case: If the biases of the samples are free pa- 
rameters to marginalize over we find that constrains on 7 degrade 
only slightly compared with the bias fixed case. In particular for 
the thinner redshift bins configurations. This is because biases are 
determined with relative errors smaller than 1% (bottom panel). 



the area) and for pre-determined galaxy samples (e.g LRGs, 
CMASS). In a photometric survey these aspects change rad- 
ically because there is no pre-selection (beyond some flux 
limits) and the number of sampled galaxies is typically very 
large (at the expense of course of poor redshift resolution). 
Therefore is interesting to investigate if the overall density 
of the samples have any impact in our results. 

Figure[7]shows the constrain in 7 for the combination of 
two samples, one unbiased population with cr z / (I + 2) = 0.05 
and a population with b — 2 and a z /(l + z) = 0.03. We keep 
the number density for the unbiased population as n(z = 
0.9) = 1.8 x 10~ 2 h 3 Mpc 3 while we vary the number density 
of the second (typically brighter) sampl^] The solid black 
line corresponds to the case in which both populations are 
correlated (same sky) and the dashed blue line to different 
areas. In both scenarios we see that decreasing the number 

9 Note that we assume the same shape for N(z) as given in Eq. |TJ 
but we vary the overall normalization, which we characterize by 
the comoving number density at 2 = 0.9. 
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Figure 7. The impact of shot noise: We consider the combined 
analysis of two populations in a rcdshift bin configuration of 
+ z) = 0.05 and show how constrains on 7 depend on 
the (shot) noise level of the more biased population (typically 
the brighter, less abundant sample). Constrains are almost un- 
affected unless the density drops by an order of magnitude or 
more compared to the one of Popl (n 2 = 0.023 h 3 Mpc" 3 ). 



density of the second population does not impact the error 
on 7 unless one degrades it by an order of magnitude or 
more (below n(z = 0.9) ~ 3.0 • 10~ J ). Above this value, the 
error is mostly controlled by the tracer with lower bias. 



3.3 Constraining the redshift evolution of the 
Growth Rate of Structure 

So far we have used the combined analysis of all the redshift 
bins to constrain one global parameter, namely the growth 
rate index 7 in Eq. (151. We now turn into constraining 
f(z)as(z) itself, as a function of redshift. We use a redshift 
bin configuration given by Az/(1 + z) = 0.1, in the photo- 
metric range 0.4 < z < 1.4. This configuration consist of 6 
bins, and hence we fit f(z)as(z) evaluated at the mean of 
these bins. These fas values are of course correlated, and 
we include the proper covariance among the measurements 
(i.e. we do a global fit to the 6 values simultaneously). 

In the left panel of Fig. [8] we focus on the gain from 
adding cross-correlations among the bins, and show the con- 
strain on fas for a single unbiased population with photo- 
metric redshift of a z = 0.05 (Pop 1, in blue) and also for a 
single tracer with bias 6 = 2 and a z = 0.03 (Pop 2, in red). 
Dashed lines corresponds to using only auto-correlations and 
solid to including also all the redshift bins cross-correlations 
to the observables. The trend for the errors when we only 
use auto-correlations are similar to the ones observed in Fig. 
8 of Ross et al. 2011. Although in detail we are using differ- 
ent widths for our redshift bins, and we use Ce while they 
used angular correlation functions, w(0). 

As in Sec. |3.1| there is a gain from the addition of cross- 
correlations, which is now split across the bins (i.e. 20 — 30% 
for Popl in each of the 6 bins, and a bit less for Pop2). 

In turn, the right panel of Fig. [8] focuses in the gain 
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Table 2. Error in the growth rate 7 from a combination of 21 
narrow bins in the range 0.94 < 2 < 1.06. The 4 top entries 
correspond to a Survey with Broad-Band (BB) filters: Popl-BB 
assumes b = 1 and <t z /(1 + 2) = 0.05 ("main sample") while 
Pop2-BB has b = 2 and a z /(l + 2) = 0.03 ("LRG sample"). The 
4 bottom entries correspond to a Survey with Narrow-Band (NB) 
filters. Here Popl and Pop2 have the same bias as the BB case 
but much precise photo-z, both with cr z /(l + 2) = 0.003. 



from combining the two tracers (and using both auto and 



cross correlations among redshift bins, as in Sec. 3.2|. Here 
the solid lines correspond to the single population cases dis- 
cussed above, while the black short-dashed line to the com- 
bined analysis assuming these populations are correlated 
(same sky). For completeness the dashed green line is the 
result when these two samples are assumed independent. 
Again, there is a factor of ~ 2.5 to be gained by combining 
galaxy samples as opposed to only the unbiased sample. 

If we compare our predictions to measurements from 
spectroscopic surveys like VIPERS (de la Torre et al. 2013) 
with constrains fa 8 (z = 0.8) = 0.47 ± 0.08 or WiggleZ 
(Blake et al. 2011) where fa s {z = 0.76) = 0.38 ± 0.04 we 
find that DES can achieve the same level of errors (~ 15%) 
in determining the growth of structure but extending the 
constrains beyond redshift of unity. This is quite unique and 
interesting as there is, to our knowledge, no other spectro- 
scopic survey expected to provide such measurements in the 
medium term future (before ESA/Euclid or BigBOSS). 

As we did for 7, we checked that when we marginalize 
over photometric errors we find differences smaller than 1% 
in the recovered constrains in f(z)a&(z) with respect to the 
case in which we assume perfect knowledge of the redshift 
distributions. For concreteness we did this cross-check for 
the case Popl x Pop2 in last two redshift bins shown in Fig. [8] 



3.4 The case of high-photometric accuracy 

In the previous sections we have focused in galaxy surveys 
with broad-band photometry for which the typical photo- 
metric error achieved is of the order 0.1 depending on galaxy 
sample and redshift (we assumed 0.03 — 0.05 (1+z)). We now 
turn to narrow-band photometric surveys such as the ongo- 
ing PAU or J-PAS Surveys (Benitez et al. 2009; Gaztanaga 
et al. 2012; Taylor et al. 2013) for which the typical ra- 
dial accuracy is a factor of 10 times better: ~ 0.003(1 + z) 
(or 10/i _1 Mpc). This scenario then resambles quite closely 
a purely spectroscopic survey (Asorey et al. 2012). 

We again study two populations, one corresponding to 
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Figure 8. Constrains on fag: derived at different redshift bins, for a bin configuration of Az/(1 + Z) = 0.1. The left panel focuses 
on one population only fits, and the gain from using auto+cross correlations among all redshift bins as observables instead of just the 
auto-correlations. The right panel stresses instead the gain from combining two populations (through their auto and cross-correlations) 
either in different patches of the sky (Popl+Pop2) or the same (PoplxPop2). In all cases Popl refers to a galaxy population with 6 = 1 
and <t z /(1 z) 0.05 and Pop2 to b = 2 and <r z /(l + z) = 0.03. The covariance among the derived errors on f(z) X 0-8(2:) is taken into 
account in the fit. Our results show that by using RSD with two tracers a DES-like photometric survey can place ~ 15% constrains in 
the evolution of fas for several bins in z > 0.8 (with errors almost uncorrelated between bins, see text for details). 



the main sample with bias 6=1 and another to the LRG 
sample with 6 = 2, both with a very good photometric ac- 
curacy of cr z /(l + 2) = 0.003. We consider a set of 21 narrow 
redshift bins of width Az = 0.003(1 + z) concentrated in 
0.94 < z < 1.06 (hence we are only looking at a portion of 
the survey redshift range). 

The error on 7 are given in Table [2] for both the new 
narrow-band and the broad-band samples discussed previ- 
ously. For a single population, this table shows that a factor 
of ~ 10 better a z yields a factor of ~ 10 gain in constraining 
power. The improvent in 7 seems to increase linear with the 
improvement in a z . 

After combining the two populations we see that the er- 
rors in 7 for the broad-band case is similar if samples cover 
the same region of sky (Popl x Pop2) or different regions 
(Popl + Pop2). This is because the redshift range considered 
(0.94 < 2 < 1.06) is too narrow compared to <r z and the cos- 
mic variance cancelation can not take place. Instead, for the 
narrow band surveys we find a 43% improvement for the case 
Popl x Pop2 with respect to Popl + Pop2. For the same sky 
case, the final error is A7 ~ 0.0163 x (5000 deg 2 /Area) 1/2 , 
in such a way that even a moderate survey of 250 deg 2 could 
achieve A7 ~ 0.07. In that same narrow redshift range, DES 
yields an error 5 times worse with 20 times better area (but 
note that in the case of small areas we could be limited by 
the i m i n , the largest scales available). 



4 CONCLUSIONS 

We have studied how measurement of redshift-space distor- 
tions (RSD) in wide field photometric surveys produce con- 
strains on the growth of structure, in the linear regime. We 
focused in survey specifications similar to those of the on- 
going DES or PanSTARRS, that is, covering about 1/8 of 
sky up to z ~ 1.4, and targeting galaxy samples with photo- 
metric redshift accuracies of 0.03 — 0.05(1 + z) (and hundred 
of million galaxies prior to sample selection). We also show 
results for ongoing photometric surveys, such as PAU and 
J-PAS, that have a much better photometric accuracy. 

First, we have found that for a single population we 
can reduce the errors in half by including all the cross- 
correlations between radial shells in the analysis. This is 
because one includes large scale radial information that was 
missed when only considering the auto-correlations of each 
bin. The final constraining power depends on the details of 
the population under consideration, in particular the bias 
and the photometric accuracy. Less bias gives more relative 
importance to RSD in the clustering amplitudes. In turn, 
better photo-z allows for narrower binning in the analysis 
and more radial information. We find that the 7 constrains 
depend roughly linearly in both bias or a z . This means that 
for 10 times better photo-z errors, such as in PAU, we can 
improve by 10 the cosmological constrains. 

Typically less bias implies a fainter sample, with worse 
photo-z, therefore these quantities compite in determin- 
ing the optimal sample. Furthermore we find that opti- 
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mal constrains are achieved for bin configurations such that 
Az ~ a z . Although the optimal errors depend on the details 
of the galaxy sample and binning strategy, the gains from 
adding cross-correlations are very robust in front of these 
variations. 

In order to avoid sample variance, we have also consid- 
ered what happens if we combine the measurement of RSD 
using two different tracers. This is motivated by the idea put 
forward in McDonald & Seljak 2009 for the case of spectro- 
scopic (hence 3D) redshift surveys, where the over-sampling 
of (radial + transverse) modes allows a much better preci- 
sion in growth rate constrains, as long as samples are in the 
low shot-noise limit. Combining auto and cross angular cor- 
relations in redshift bins, we find that if we assume that both 
tracers are independent, which corresponds to samples from 
different regions on the sky, the constrains on the growth 
of structure parameters improve a 30-50% (due to the fact 
that one has doubled the area). Remarkably if we consider 
that the populations are not independent, i.e., they trace the 
same field region, we find an overall improvement of ~ 2 — 3 
with respect to single populations when constraining 7. This 
means that there is a large potential gain when sampling the 
same modes more than once. 

Translating into actual constrains this implies that a 
DES-like photometric survey should be able to measure the 
growth rate of structure 7 to an accuracy of 5 — 10% from 
the combination of two populations and all the auto+cross 
correlations in the range 0.4 < z < 1.4 (see Fig. [T|. Even 
though these values correspond to a survey of 5000deg 2 
(fsky = 0.125) they should scale as f~^J 2 for a different 



area, given our assumptions for the covariance in Eq. ( 12 \ 



In Fig. [7] we have shown that constrains weaken once 
one of the populations enter a shot-noise dominated regime, 
as is typical of spectroscopic samples. However one needs to 
dilute over 10 times the number densities for a photometric 
survey, such as DES, for this to happen. Thus, as shown in 
Section 3.4, by improving on photo-z accuracy without much 
lost of completeness, a photometric sample can in fact out- 
perform a diluted spectroscopic version with similar depth 
and area (see also Gaztanaga et al. 2012). In this paper 
we focused on large angular scales where the approximation 
of linear and deterministic bias and linear RSD should hold 
(see for instance Crocce, Cabre, & Gaztanaga 2011). Al- 
though we set tmax ~ 200, much of the constraining power 
in our results, given the typical size of our redshift bins, 
comes from larger scales, I < 40. Yet, a more realistic as- 
sessment of these aspects will need to resort to numerical 
simulations. We leave this for future work. 

Lastly, we also investigated what constrains can be 
placed with this method in the evolution of the growth rate 
of structure, f(z) x crg(z). We found that binning two DES 
populations into 6 bins in the range 0.4 < z < 1.4 yields 
constrains on f(z) x as(z) of ~ 15% for each bin above 
z ~ 0.6. That case corresponded to bin widths larger than 
the photometric errors of the samples, which may not be op- 
timal but yield constrains almost uncorrelated between bins 
(Pij ~ — 0.05^p"] Figure [9] shows instead the results from a 

10 For bins we define the cross-correlation coefficient pij as 
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Figure 9. Combined constrains in the evolution of the growth 
rate of structure from spectroscopic data, 2dFGRS, SDSS-LRGs, 
Wiggle-Z and BOSS (see text for details), and forecasted for DES 
using two photometric populations (same as in Fig. 8). The ad- 
dition of DES (shaded area) allows to trace the growth rate of 
structure all the way to z ~ 1.4. 



narrower binning, Az/(1 + z) — 0.05. This leads to better 
constrains, A (/erg) ~ 10%, at the expense of more correla- 
tion between bins, 0.2 < pij < 0.65. 

In addition to the DES forecast (shadowed region) we 
over-plot in Fig. [9] current constrain from spectroscopic sur- 
veys, 2dFGRS (Percival et al. 2004), LRG's from SDSS 
(Tegmark 2006 and Cabre & Gaztanaga 2009), WiggleZ 
either from power spectrum (Blake et al. 2011) or correla- 
tion function (Contreras et al. 2013) and the recent BOSS 
results (Reid et al. 2012). Note that these values are not ex- 
pected to improve radically in the near future. This implies 
that DES will be able to add quite competitive constrains 
in a redshift regime unexplored otherwise with spectroscopic 
surveys (i.e. z > 1), yielding a valuable redshift leverage for 
undertanding the nature of dark energy and cosmic acceler- 
ation through the growth of structure. 
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